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•n . Whereas ground state spectroscopy for quenched QCD is well understood, it is still a chal- 
^ ! lenge to obtain results for excited hadron states. In our study we present results from a new ap- 
proach for determining spatially optimized operators for lattice spectroscopy of excited hadrons. 
In order to be able to approach physical quark masses we work with the chirally improved Dirac 
O ' operator, i.e., approximate Ginsparg-Wilson fermions. Since these are computationally expen- 
sive we restrict ourselves to a few quark sources. We use Jacobi smeared quark sources with 
different widths and combine them to construct hadron operators with different spatial wave 
I functions. This allows us to identify the Roper state and other excited baryons, also in the 
^ I strange sector. 

^ 1. Quark sources and interpolating fields 

c3 ■ Ground state spectroscopy on the lattice is by now a well understood physical problem with 
impressive agreement with experiment. The lattice study of excited states is not as advanced. 
The Euclidean correlation function of an interpolating operator contains contributions from all 
states with the correct quantum numbers and the masses of the excited states show up only in 
the sub-leading exponentials. A direct fit to a sum of exponentials is cumbersome since the 
signal is strongly dominated by the ground state. Also with methods such as constrained fits, 
the maximum entropy method and a recently proposed method one still needs high 

statistics for reliable results. 

An alternative method is the variational method [ ^\ where one diagonalizes a matrix con- 
taining all cross-correlations of a set of several basis operators Oi with the correct quantum 
numbers, 

C(t)., = {0,{t) 0]m = Y: (0|O»e-**^"(n|O]|0) , (1) 

n 

(written here for infinite temporal extent of the lattice and operators projected to vanishing spa- 
tial momentum). For a large enough and properly chosen set of basis operators each eigenmode 
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is then dominated by a different physical state. Finding these states is equivalent to solving the 
generalized eigenvalue problem [,4J with eigenvalues behaving as 

X^'\t) oc e-'^'^ [l + 0(e"*^*^^)] . (2) 

Each eigenvalue corresponds to a different energy level dominating its exponential decay; 
AMfc is the distance to nearby energy levels. The optimal operators which have maximal over- 
lap with the physical states are linear combinations of the basis operators. After normalization 
the largest eigenvalue gives the correlator of the ground state, the second-largest eigenvalue 
corresponds to the first excited state, and so on. 

Although the set of basis operators can be made arbitrarily large, in realistic calculations the 
intrinsic statistical errors make the diagonalization increasingly unstable. The challenge is then 
to provide a set of operators large enough to span the physically relevant space but still small 
enough to allow significant results for the given statistics. 

In earlier work [ 5 1 we used three different interpolating fields for the nucleon sector, but 
could not resolve a signal for the Roper excitation. Here we concentrate on the first two of 
those, 

Xi{x) = eabc [u1{x)C-f5db{x)\ Uc{x) , X2ix) = eabc [u1{x)Cdb{x)\ j^Ucix) . (3) 

We concluded that in order to get a reliable Roper signal it is also important to optimize 
the spatial properties of the interpolating operators. It can be argued that a node in the radial 
wave function is necessary to capture reliably the Roper state or other radially excited hadrons. 
In [ we demonstrated that an elegant solution is to combine Jacobi smeared quark sources 
with different widths to build the hadron operators and compute the cross-correlations in the 
variational method. We find good effective mass plateaus for the first and partly the second 
and even the third radially excited states. The eigenvalues X'^^\t) are then fitted using standard 
techniques. 

Already earlier [ |7| Jacobi smeared sources were combined with point sources and cross- 
correlations studied in similar spirit (see also ['HI)- The technique of Jacobi smearing is well 
known [|71|9|. The smeared source lives in the timeslice t = and is constructed by iterated 
multiplication with a smearing operator H ona point-like source. The operator H is the spatial 
hopping part of the Wilson term at timeslice 0; it is trivial in Dirac space and acts only on the 
color indices. This construction has two free parameters: The number of smearing steps and 
the hopping parameter k. These can be used to adjust the profile of the source. We worked with 
two different sources, a narrow source n and a wide source w with parameters chosen such that 
the profiles approximate Gaussian distributions [|6|. The distribution of the two source-types 
correspond to half-widths of a ~ 0.27 fm and a ~ 0.41 fm. 

The two sources allow the system to build up radial wave functions with and without a node. 
The parameters were chosen such that simple linear combinations c„ n -|- c^jW of the narrow 
and wide profile can approximate the first and second radial wave functions of the spherical 
harmonic oscillator: The coefficients c„ ~ 0.6, ~ 0.4 approximate a Gaussian with a half- 
width of cr ~ 0.33 fm, while c„ ~ 2.2, ~ —1.2 approximate the corresponding excited wave 
function with one node. Note, however, that the final contribution of the corresponding source 
terms is not pre-determined but results from the diagonalization. 

Following the variational method [ 4 1 we compute the complete correlation matrix of opera- 
tors. The hadron sources we use for the correlation matrix are constructed from the narrow and 
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Figure 1. Masses of S and S baryons in units of the lattice spacing a ~ 0.148 fm (=1/1.33 GeV) 
(triangles: ground state, circles: negative parity ground state, squares: first positive parity ex- 
cited state, full symbol: experimental value of ground state); the mass of the strange quark has 
been adjusted to give physical iC-meson mass. 



wide quark sources. Respecting symmetries we then have for both types of nucleon operators 
dSj six combinations of wide or narrow quark sources, in total 12 interpolating field operators 
and a 12 X 12 correlation matrix. The final form of the wave function is determined through the 
eigenvectors resulting from the variational method. 



2. Baryon results 

We use the chirally improved Dirac operator [[TOl| and work in the quenched limit, ie., without 
dynamical fermions. This operator is an approximate Ginsparg- Wilson operator and was well 
tested in quenched ground state spectroscopy [|TT1 where pion masses down to 250 MeV have 
been reached at a considerably smaller numerical cost than needed for exact Ginsparg-Wilson 
fermions. The gauge configurations were generated on a 12^ x 24 lattice with the Liischer-Weisz 
action [ 12 1 (further studies for larger lattices are in progress). The inverse gauge coupling is 
P = 7.9, giving rise to a lattice spacing of a = 0.148(2) fm as determined from the Sommer 
parameter [\T3\. The statistics of our ensemble is 100 configurations. For the u, d quarks we 
use degenerate quark masses m ranging from am = 0.02 to am = 0.20. More details are 
discussed in [ 6|. 

Although we determine the full correlation matrix as discussed, it turns out to be numerically 
advantageous to analyze only a subset of operators. The exponential decay of the three leading 
eigenvalues is clearly identified. We identify these signals with the nucleon, the Roper state and 
the next positive parity resonance A^(1710). The larger eigenvalues have too large statistical 
errors for a reliable interpretation (see also [ OJ concerning the problem of nucleon-?]' ghost 
contributions). 

In principle both interpolators xi and X2 should couple to all states with the corresponding 
quantum numbers; however, the coupling amplitude depends on the internal structure of the 
physical states (with possibly different overlap). Since we can continuously increase the quark 
mass and thus make contact with the heavy quark region we also may argue in such a framework. 
All Roper states form an excited 56-plet of SU(6) and in this multiplet the parity of any two- 



4 



quark subsystem is positive. The A^(1710) belongs to another muMplet containing both positive 
and negative parity two-quark subsystems. The two-quark subsystem in the brackets in xi has 
positive parity, while the subsystem in the brackets of X2 has negative parity. Hence in X2 
we should see only the signal from the A^(1710) state, and no signals from the nucleon and 
the Roper. This is exactly what we observe in [|6|, where the relevant figures for the nucleon 
channel are shown. 

We may also combine quarks with different masses, allowing for a strange quark. We fix the 
strange quark mass such that the K-meson mass has its experimental value in the chiral limit 
of the u, d-masses. We then can determine also the S- and S-baryons using interpolator xi and 
again different width quarks sources. First results are exhibited in Fig. [T] We have to stress, 
however, that the results given here are for relatively small lattices (spatial lattice size 1.8 fm). 

A crucial test of our method is to check whether indeed the ground state is built from a 
nodeless combination of our sources and the excited states do show nodes. This question can 
be addressed by analyzing the eigenvectors of the correlation matrix. This has been done in 
and indeed confirms the expectation. 

3. Conclusion 

Ground state masses approach, in naive chiral extrapolation, their experimental values well. 
The excited state masses are above their experimental values. Since our results are for small 
spatial lattice size a volume dependence study still has to confirm this behavior. Excited states 
may be more affected by finite volume effects, but also be more sensitive to the quenched 
approximation. Finite volume studies, including also scaling properties, are under way. 
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